Improved Empirical Parametrization of Fragmentation Cross Sections 
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A new version is proposed for the universal empirical formula, EPAX, which describes fragmen- 
tation cross sections in high-energy heavy-ion reactions. The new version, EPAX 3, is shown to 
yield cross sections that are in better agreement with experimental data for the most neutron-rich 
fragments than the previous version. At the same time, the very good agreement of EPAX 2 with 
data on the neutron-deficient side has been largely maintained. Comparison with measured cross 
sections show that the bulk of the data is reproduced within a factor of about 2, for cross sections 
down to the pico-barn range. 

PACS numbers: PACS: 25.70.Mn, 25.70.Pq, 25.75.-q 



INTRODUCTION 



During the last two decades, projectile fragmentation 
and separation at high energies has become one of the 
most important methods for the production of exotic nu- 
clei. Production and separation of such species can be 
performed at energies of around 500 to 1000 A MeV 
at the SIS/FRS facility Q at GSI in Darmstadt, Ger- 
many. At PJKEN in Wako, Japan, the RIBF facility 
with the BigPJPS separator [2] allows to perform such 
experiments at energies of up to 350 A MeV. Consider- 
ably lower energies of up to 140 A MeV are attained at 
the A1900 separator Q at NSCL in East Lansing, USA, 
whereas the GANIL/LISE facility at Caen, France, op- 
erates below about 80 A MeV In all cases, the main 
advantage of high-energy projectile fragmentation as a 
production method lies in the fact that very clean frag- 
ment beams can be produced; the contaminants become 
less and less important with increasing energy. Moreover, 
the high energies allow unambiguous event-by-event iso- 
tope identification of the fragments produced. Again, 
the higher energies help in suppressing unwanted charge 
states of the fragments. 

All of the above-mentioned facilities rely on good es- 
timates of fragmentation cross sections in order to pre- 
dict the production rates of exotic nuclei, together with 
ion-optical transmission calculations with codes like MO- 
CADI 0] or LISE++ @| . One important tool to predict 
fragmentation cross sections is a universal analytical for- 
mula called EPAX @, i|. This formula allows to calcu- 
late the yields from fragmenting all non-fissile projectiles 
in the range of projectile masses, A p , between about 40 
to 209. In particular, the formula tries to take prop- 
erly into account the influence of the projectile proton 
or neutron excess onto the neutron-to-proton ratio of the 
fragments. The cross sections are assumed to be energy- 
independent, which seems to be supported by most of the 
measured data, even though no systematic studies of the 
bombarding-energy dependence have been performed for 
targets like 9 Be. The EPAX parametrization aims at re- 
producing the bulk of the measured cross sections within 
a factor of two for fragment masses down to about half 



the projectile mass. This latter restriction is no severe 
limitation since the highest transmission in an ion-optical 
projectile- fragment separator is always attained for small 
mass and charge differences between projectiles and frag- 
ments. 

The previous version of the formula, EPAX 2 Q, has 
proven to give rather realistic estimates of many pro- 
duction cross sections of exotic nuclei. As an example, 
the very proton- rich fragment 100 Sn was produced with 
11.2 ±4.7 pb @ in 124 Xe+ 9 Be (later remeasured to yield 
5.8±2.1 pb by Straub et al. 10]), not far from the EPAX 
2 prediction of 7.4 pb. On the other hand, large dis- 
crepancies with measured data were found for extremely 
neutron- rich nuclei like, e.g., fragment yields from 1 A 
GeV 136 Xe+ 9 Be, which were overestimated by EPAX by 
up to two orders of magnitude (Tl| . It is therefore desir- 
able to try and modify the EPAX 2 parameters in such 
a way that also the yields of very neutron-rich fragments 
of medium- and heavy-mass projectiles can be predicted 
with better accuracy. 

In the following, the basic characteristics of the EPAX 
formula will be reviewed. This is followed by an expla- 
nation how the parameters of EPAX can be obtained 
by fits to recent data sets. These fits allow to get a 
feeling how the EPAX parameters vary with fragment 
mass and with proton or neutron excess of the projectile. 
The final choice of the parameters for EPAX 3 is then 
presented, together with comparisons between measured 
data of selected systems and the predictions of the new 
parametrization. 



II. THE EPAX FORMULA 

As explained in detail in Ref. Q, the cross section a 
of a fragment with mass A and charge Z produced by 
projectile fragmentation from a projectile {A p ,Z p ) im- 
pinging on a target {At, Z t ) is written as 

a(A, Z) = Y A crz(Zprob-Z) = Y A n exp (-R \Z prob -Z\ u ) 

The first term, Ya, represents the mass yield, i.e. the 
sum of all isobaric cross sections with fragment mass A. 
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The second term, az, describes the "charge dispersion", 
the distribution of elemental cross sections with a given 
mass, A, around its maximum, Z pro i,. The shape of the 
charge dispersion is controlled by the width parameter, 
R, and the exponent, U. The factor n = U R/ir simply 
serves to normalize the integral of the charge dispersion 
approximately to unity. Note that the isobar distribu- 
tions are not symmetric on the neutron- and proton-rich 
side, therefore U has two different values, U p and U n , 
on the proton- and neutron-rich side of the valley of /in- 
stability, respectively. In Ref. Q , the exponent U for the 
neutron-rich side of the isobar distribution was chosen as 
U n = 1.65. In the present version, the exponent on the 
proton-rich side is also taken as a constant, U p = 2.1, 
close to the value of 2.0 for a Gaussian. Thus, in a 
fitting procedure of the EPAX parameters to measured 
data, three parameters have to be obtained for each frag- 
ment mass, A: An amplitude constant, proportional to 
the mass yield, Ya, the centroid of the charge disper- 
sion, Zprgb, and the width parameter, R. As was shown 
in Ref. [8], a fourth parameter has to be determined for 
those data sets where small cross sections of very proton- 
rich fragments have been measured: the transition point, 
Z e xp, where the quasi-Gaussian shape of the charge dis- 
persion on the proton-rich side turns into an exponential 
decay. 



A. EPAX fits to measured data sets 

During the decade between the release of EPAX 2 
and today, many data-sets of experimental fragmentation 
cross-sections have been published, some of them very 
comprehensive in their coverage of fragment charge and 
mass numbers. In addition to the experiments cited in 
the EPAX 2 paper [§[ and the ones mentioned above 0- 
[Tl| . very comprehensive data on the fragmentation of 
40 > 48 Ca and 58 > 64 Ni on 9 Be and 181 Ta targets at 140 A 
MeV by Mocko et al. [HJ turned out very useful. Tarasov 
et al. 13] have published a very careful study of 76 Ge 
fragmentation at 132 A MeV, where great care was taken 
to reliably model the (relatively small) transmission of 
the MSU A1900 separator at this bombarding energy. 
Data for heavier projectiles, all on 9 Be targets, were mea- 
sured at the FRS separator at GSI Q: For 92 Mo (at 500 
A MeV) by Fernandez et al. [3, for 112 Sn by Stolz et 
al. 0, and for 208 Pb by Benlliure et al. 0, the latter 
all measured at 1 A GeV. In addition, a very compre- 
hensive data set for 124,136 Xe on a 208 Pb target, again 
at 1 A GeV, has been published by Henzlova et al. [18j . 
A remarkable milestone was reached by Perez and co- 
workers [l|| who succeeded to measure a rather large set 
of fragmentation cross sections for a secondary beam of 
the extremely neutron- rich nucleus 132 Sn, obtained from 
238 U projectile fragmentation at 1 A GeV. 

In a few cases, cross sections for the same fragment 
have been measured at both, medium and high bom- 
barding energies. A very comprehensive data set is avail- 



able for the projectile 58 Ni, with data measured at 140 
A MeV [3 and at 650 A MeV The experimen- 

tal cross sections agree well within their respective error 
bars, therefore the 140 A MeV data sets for 40,48 Ca and 
58,64 Ni on 9 Be of Ref _ rjjj wnere i nc l u ded in the fits of 
Eq. (fTJ) to isobaric cross section distributions. 

The EPAX formula, Eq. ([1]), assumes that the free pa- 
rameters of the formula, Ya, Z pro i,, and R, vary smoothly 
with fragment mass, A (the parameter U has a fixed value 
of Up = 2.1 on the proton-rich side and of U n = 1.65 
on the neutron-rich side). To obtain a feeling for the 
A-dependence of the parameters, least-squares fits of 
Eq. (TT]) with 3 parameters were made to the cross sec- 
tions of many A chains of the various systems mentioned 
above. Only for 58 Ni, the smallest cross sections on the 
proton-rich side required the inclusion of a fourth pa- 
rameter, Z exp , which will be explained in detail in Sub- 
sect. Ill B 41 It turned out that Eq. (fTJ) can describe the 
distributions with very good accuracy over many orders 
of magnitude. The significance of a good fit is high, how- 
ever, only for distributions where many data points on 
both, the proton and neutron rich side are available. If 
only few data points are available for a given A and, in 
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FIG. 1. (color online) Experimental cross sections for A = 50 
isobars from 58 Ni fragmentation at 140 A MeV (Ref. 0, 
circles) and at 650 A MeV (Ref. [2(j, squares) in a 9 Be target. 
The full curve represents a least-squares fit of Eq. {1} to the 
data. The dashed curve shows the agreement of the EPAX 3 
formula with the data. 

particular, if only one side of the charge dispersion has 
been measured, it is difficult to obtain a unique set of fit 
parameters since the parameters are strongly correlated. 

A good example is visualized in Fig. [TJ where Eq. (fTJ) 
has been fitted to data for 9 different isobars of the 
A = 50 mass chain, where 3 cross sections have been 
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measured at both energies, 140 A MeV [l2| and at 650 A 
MeV (20| . The smallest cross sections measured amount 
to sub-nanobarns. Note that on the proton-rich side the 
quasi-Gaussian with an exponent of U p — 2.1 turns to an 
exponential slope at Z exp — 25.04 (see Subscct. Ill B 41) . 

Fits of this type were made for all of the systems of 
Refs. [13, [H M, M HI for those isobaric chains 
where at least 4 data points are available. For 208 Pb, 
the data measured with a nat Cu target by De Jong et 
al. [2l| were combined with those measured with a 9 Be 
target by Benlliure et al. [l6[ by scaling the former ones 
with a scaling factor of 0.77 to account for the larger size 
of the Cu target (see Eq. (|3|) below). Examples of the 
fragment- mass dependence of the parameters Y A , Z pro b, 
and R are given below in Figs. [Mil 
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B. The EPAX 3 formula and its parameters 

1. Mass yields 

As noted already in Rcf. [Jj, the mass-yield curve, 
Y A (A), can to a large extent be approximated by an ex- 
ponential depending on the difference between projectile 
and fragment masses, (A p — A). The slope of this expo- 
nential, P, is a function of the projectile mass only. The 
functional form of the mass yield curve in EPAX 3 has 
been kept identical to the previous version: 

Y A = S-P-exp(-P-(A p -A)), (2) 
S = Sl ■ (A 1 / 3 + A] /3 - s 2 ) [barn], (3) 
P = exp(pi+p 2 -A p ). (4) 

Eq. (j3|) shows that the dependence of S on the sum of 

1 /3 1 /3 

A t and A t and thus on the circumference of the col- 
liding nuclei has been kept. As the only change the over- 
lap factor, S2, has been slightly reduced from 2.38 to 
1.8 fm. Results from a recent paper on charge-changing 
cross sections induced by 650 A MeV 40 Ar and 1 A GeV 
48 Ti ions by Zeitlin et al. [l7| confirm that Eq. © is a 
good description of the geometrical scaling factor. Fig. [5] 
shows the ratios of charge-changing cross sections for Al, 
Cu, Sn, and Pb targets, relative to those of a carbon tar- 
get, and averaged over all charge changes measured in 
Ref. [I?]]- The x-axis corresponds to the scaling variable 
in Eq. ©, (Ap^ 3 + A]^ 3 — s 2 ); the straight line represents 
Eq. ©. 

In Ref. 8] it was shown that the mass yield Y A has to 
be increased near the projectile (see Eq. (13) of Ref. Q). 
The same formula was kept for the new Version 3: 

Y' A = Y A - 1^ ■ A p ■ (x - yi ) 2 } (5) 

for x > yi ; x is the ratio of fragment mass to projectile 
mass, x = A/Ap. The numerical values of all constants 
can be found in Table |U 

Typical examples of mass yields are shown in Fig. [3] 
The data points represent the numerical values of Y A as 



FIG. 2. Charge-changing cross sections of relativistic Ar 
and 48 Ti beams interacting with the targets indicated in the 
figure, relative to those with C targets. The ratios were av- 
eraged over all charge changes measured in Ref. JT3]. The 
straight line represents the EPAX description, Eq. ©. 

obtained from fitting Eq. (Q} to isobaric chains of cross 
sections of Refs. [12|, |l8|. One can clearly see that the 
slopes, P, are much steeper for lighter projectiles like 
58,64 Ni than for heavier ones like 124,136 Xe, corroborat- 
ing the A p dependence of P, Eq. Eq. ([5]) has been 
chosen such that the upward bend of Y A towards the 
projectile mass scales with A p and thus has only a minor 
effect on the 58 > 64 Ni mass yield distributions. The slope 
differences between each of the two Ni and Xe isotope 
pairs are entirely due to the ^-dependence of P. 
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FIG. 3. (color online) Mass yields, Y A (A), for 58 ' 64 Ni+ 9 Be 
and 124 ' 136 Xe+ 208 Pb. The open symbols represent the nu- 
merical values of Y A derived from fits of Eq. |T|) to the ex- 
perimental data of Refs. p3, [3- The full lines have been 
calculated from Eqs. ([2][5]) 
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2. Centroid Z pro b 

As mentioned above, the charge dispersion is mainly 
characterized by its centroid, Z pro b, and its width pa- 
rameter, R. As in EPAX Versions 1 and 2, the loci of 
maximum cross section, Z pro b(A), have been parameter- 
ized relative to the valley of /3-stability, 

Zprob = Zp + A + A m . (6) 

Zp is approximated by the smooth function 

Zp = ,4/(1.98 + 0.0155 -A 2/3 ). (7) 

A is found to be a linear function of the fragment mass, 
A, for heavy fragments (A > di), and is extrapolated 
quadratically to zero: 

(ch + dz-A ifA>d 4 

\d 3 - A 2 if A < d A y ' 

The sum Za + A defines the "residue corridor" , the locus 
in the A — Z plane where the fragments from all projec- 
tiles located near /3-stability end up after long evapora- 
tion chains. This corridor is close to the line Charity [22[ 
has found from studying the residue distributions after 
evaporation from highly excited compound nuclei (see 



TABLE I. Constants used in the new EPAX 3 formula. 



Parameter Constant Value 

Scaling factor S si 0.27 

s 2 1.80 

Mass yield slope P pi —1.731 

p 2 -0.01399 

Mass yield correction factor j/i 0.75 

y 2 0.10 

Z prob shift A di -1.087 

d 2 3.047 • 10" 2 

d-j, 2.135 • 10~ 4 

d 4 71.35 

d 5 -25.0 

d 6 0.80 

n-rich memory effect m 0.40 

n 2 0.60 

p-rich memory effect AJ^ qi —10.25 

q 2 10.25 

n-rich slope U n U n 1.65 

p-rich slope U p U p 2.10 

Width parameter R ro 2.78 

ri -0.015 

r 2 3.20 ■ 10~ 5 

r 3 0.0412 

r 4 0.124 

r 5 30.0 

r 6 0.85 

Transition point to h 1.20 

exponential slope l 2 0.647 

Downscale factor b\ 2.3 • 10 -3 

b 2 2.4 




fragment mass A 



FIG. 4. (color online) Loci of the peak cross sections, 
Zp ro b(A), as a function of fragment mass, for various projec- 
tiles as indicated in the Figure. For clarity, the value of Zp(A) 
has been subtracted at each mass number A. The data points 
represent numerical values of (Zp — Z pro b) obtained by fitting 
Eq. {TJ to isobaric distributions for the systems indicated. 
Positive (negative) ordinate values correspond to proton-rich 
(neutron-rich) nuclei relative to Zp. The large dots indicate 
where the respective projectile is located in this graph. The 
lines represent the EPAX predictions according to Eqs. (6-11). 

dashed line in Fig. 4 and Eqs. (6,7) in Ref. [H|]). This 
points to the fact that the fragment distributions far away 
from the projectiles are largely controlled by evaporation 
from equilibrated pre-fragments. 

It is obvious that close to the projectile, Z pro b has to 
move away from the corridor and approach Z pi i.e., for 
projectiles located on the line of /3-stability, A has to 
approach zero. That has been achieved by the following 
modification 

A' = A-[l + d 5 -(x-d 6 ) 2 } (9) 

if x > dg. 

For projectiles on the proton- or neutron-rich side of 
the valley of /3-stability, an additional correction of Z pro b 
is necessary, which accounts for the "memory effect" , the 
gradually vanishing proton or neutron excess of the frag- 
ment reflecting the proton or neutron excess of the pro- 
jectile, (Z p — Zp p ), where Zp p is the value of Eq. Q 
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for the projectile, A = A p . This additional correction, 
termed A m in Eq. ©, should reflect the full excess at 
the projectile mass and should vanish gradually with in- 
creasing mass loss from the projectile. The functional 
forms have been kept as in the previous version, 

A™ = (m • x 2 + n 2 • x 4 ) • (Z p - Z Pp ) (10) 

for neutron-rich projectiles, and 

A ™ = ex P (<ii + <i2 ■ x) ■ {Z p - Zf jp ) (11) 

for proton-rich projectiles. In both cases, the choice of 
the parameters guarantees the limiting properties men- 
tioned above. 

Fig. @] visualizes how the combined effects of A, A', 
and A m compare to numerical values of (Zp — Z pro b) 
derived from fits of Eq. (CD) to the experimental data sets 
of Refs. [HI, [l8|, [24[. The zero line corresponds to 
the line of /3-stability according to Eq. ([7]). The full lines 
have been calculated using Eqs. (fSllTTj). The agreement 
with the fitted data points is rather good close to the 
respective projectile masses; further away, one notes bad 
agreement in particular for 136 Xe. 

3. Width parameter R 

Similar to the parameter Z pro b just discussed, the 
width parameter, R, is to first order a function of frag- 
ment mass only, irrespective of the projectile. In Ref. [8[ 
it was found that the experimental i?-values can be 
approximated by a simple exponential (see Eq. (8) of 
Ref. [1]). As visualized in the following section, this yields 
too wide charge dispersions for very heavy nuclei, such 
that now a parabolic dependence in log scale has been 
chosen: 

R = Rq ' 

exp (r x ■ A + r 2 ■ A 2 ) (12) 

It was already noted in Ref. 0] that the charge disper- 
sions from neutron-rich projectiles are a bit wider than 
those from proton-rich or /^-stable ones. In the present 
version of EPAX, this has been coded in the following 
form: the parameter Rq in Eq. (|12l) depends differently 
on the proton and the neutron excess of the projectile: 

Rq = r ■ exp [r 3 • (Z p - Zp p )) (13) 
RP = r - exp [r 4 • (Z p - Zp p )} (14) 

It is again obvious that close to the projectile the widths 
of the charge dispersions should shrink since the smaller 
number of nucleons removed from the projectile also re- 
duces the variance of the mass loss. Therefore, the follow- 
ing modification has been introduced near the projectile, 
i.e. for x > r§: 

R' = R-exp[r 5 - ^fT p ■ (x - r & f]. (15) 

The Ap-dependence of the pre-factor in the exponent 
guarantees a small effect for light and a stronger one for 
heavy projectiles. 




20 40 60 80 100 120 



1 - 




60 80 100 120 140 160 180 200 



1 1 1 1 1 


64 Ni 


i i i i | i 
86 Kr 


■ i | i i i 
136 Xe ■ 


□ 


ij 


a 


A . 






it 










£'] 
A// 








Ay 


: (c) 

. i . . . 


i . . . 


i ... i . 


. . I . . . 



40 60 80 100 120 140 



fragment mass A 

FIG. 5. (color online) (a) Fragment-mass dependence of 
the width parameter, R(A), for proton-rich projectiles. The 
symbols represent numerical values of R obtained by fitting 
Eq. {TJ to the systems indicated. The dashed curves denote 
the old EPAX 2 parametrization, while the new Version 3 
(Eqs. (|12ll5p ) is represented by the full curves, (b) The same 
for projectiles near the line of /3-stability. Note the reduced 
widths (larger R values) for 208 Pb-fragments according to 
EPAX 3. (c) The same for neutron-rich projectiles. Since 
here the dashed curves are almost invisible, the dash-dotted 
line indicates the EPAX 2 curve for stable and proton-rich 
projectiles. 



Fig. [S] compares the width parameters calculated from 
Eqs. (|12H15p to numerical values of R obtained by fit- 
ting Eq. ([1} to experimental cross sections of the various 
systems. The upper part, (a), shows this comparison 
for proton-rich projectiles ( 40 Ca, 58 Ni, and 124 Xe). The 
EPAX 3 parametrization (full lines) leads to more nar- 
row charge dispersions than the EPAX 2 one (dashed 
lines). Fig. [SJd has been plotted for the projectiles 
129 Xe and 208 Pb, both located close to the line of 0- 
stability. The more narrow widths compared to the pre- 
vious parametrization [8] for 208 Pb are clearly visible. 
The lower panel, (c), displays R(A) for three neutron-rich 
projectiles, 64 Ni, 86 Kr, and 136 Xe. The parametrization 
according to EPAX 2 coincides with the data points and 
is hardly visible; there is not much difference to the Ver- 
sion 3 parametrization (full lines). The dash-dotted line 
indicates the EPAX 2 exponential for stable or proton- 
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rich projectiles; it is obvious that neutron-rich projectiles 
lead to considerably wider distributions. 



III. COMPARISONS WITH EXPERIMENTAL 
CROSS SECTIONS 



4- Modifications for very proton-rich fragments 

A major modification of EPAX 2 compared to EPAX 
f consisted in the introduction of an exponential rather 
than a quasi- Gaussian shape of the charge distribution 
beyond a certain transition point on the proton-rich 
side 0. This was coded by calculating the derivative 
of the logarithm of the cross section (Eq. (TTJ)): 

dF d(log(a)) _-2-R _ v 

^prob ) 



dz 



dZ 



(16) 



ln{10) 

The transition point to the exponential slope, Z expi has 
been calculated for the proton-rich side as a function of 
the fragment mass A according to 



Zexp(A) — Zp ro i,(A) + 



dF 
dZ 



Zra(lO) 
2 • R{A) 



(17) 



From Z exp on, the slope is exponential with the same gra- 
dient as Eq. (fT|) at this point. The empirical parametriza- 
tion for the transition slope proposed in Ref. Q, 



^- = l 1+h .(A/2r 



(18) 



not only yields perfect agreement with very proton-rich 
fragments from 58 Ni, but provides the same quality also 
for fragments from 124 Xe, e.g. for 100 Sn, for which EPAX 
2 calculates a cross section of 7.4 pb, whereas Straub et 
al. measure 5.8 ± 2.1 pb [10|. This parametrization has 
therefore been kept in EPAX 3, the 100 Sn cross section 
from this new version now amounts to 6.9 pb. 



5. Modifications for very neutron-rich fragments 

It has been noted by several authors (e.g. by Bcn- 
lliure et al. [ll], [r|) that EPAX 2 overestimates cross 
sections of very neutron rich fragments from, e.g., 136 Xe 
or 208 Pb by up to two orders of magnitude. After many 
failed attempts to modify the EPAX parameters in such a 
way that the discrepancies were removed, it became clear 
that only a "brute-force solution" could solve this prob- 
lem. This means that all cross sections from neutron-rich 
projectiles calculated by EPAX 3 with Eqs. (1-18) above 
should be downscaled by a factor that depends on both, 
the neutron excess of the projectile and the fragment. 
This "brute-force factor" can be written as 

f bf = IQiB-iZ^-Z+bif] (1Q) 

B = -h- \Z p - Z 0p \ (20) 

where Zp and Zp p are the values of Eq. ([7]) for fragment 
mass A and projectile mass A p , respectively. It applies to 
fragments where (Zp — Z) > (Z p — Zp p + b 2 ). The effect 
of this factor is quite dramatic for the most neutron-rich 
fragments of heavy projectiles as will be shown below. 



A. Projectiles close to /3-stability 

An important testing ground for the quality of EPAX 
3 predictions are the fragmentation cross sections of 
208 Pb, one of the few heavy projectiles to explore the 
yet poorly studied region in the chart of nuclides "south" 
and "south-east" of this nucleus. Data were measured at 
1 A GeV bombarding energy on nat Cu [H[ and 9 Be [HI 
targets. For the following comparison (Fig. [5]) the cross 
sections of Ref. [2l[ were scaled by a factor of 0.78 ac- 
cording to Eq. ©. Both data sets fit excellently to each 
other and are well reproduced over a large range of frag- 
ment masses by EPAX 3. In particular, the predictions 
for neutron-rich fragments farther away from the projec- 
tile are now much closer to the experimental data than 
those of EPAX 2. Small discrepancies remain for the 
most neutron-rich Tl and Hg isotopes. 



B. Fragmentation of neutron-rich projectiles 

It was mentioned in the Introduction that one of the 
most severe deficiencies of EPAX 2 was its overpredic- 
tion of cross sections of very neutron-rich fragments, e.g. 
those from 136 Xe+ 9 Be [ll| . It is therefore interesting 
to see the effect of the new parameters and in particular 
that of the "brute- force" downscale factor (Eqs. (|19I20[1 ). 
Some isotope distributions from this system are shown in 
Fig. [7] The flat isotope distribution for Z = 53 is not so 
well reproduced, but for the bulk of the cross sections of 
very neutron-rich fragments EPAX 3 fits better than the 
previous version. For the few data points measured for 
Z = 45, the new description seems even to underpredict 
slightly the data. 

Similar conclusions can be drawn if one plots the cross 
sections for the proton-loss channels with N — 82 from 
136 Xe+ 9 Be (Fig. ©. Benlliure et al. have noted in 
Ref. [ll| that EPAX 2 yields too flat a slope for the cross 
sections of the N = 82 isotones. Fig.[8]now demonstrates 
that the new parametrization gives a slope in very good 
agreement with measured data, despite a discrepancy, 
however, at Z = 53. 

Fig. [5] compares EPAX 3 also to the predictions 
of a physical model, the "cold fragmentation" model 
COFRA [25j . This model, a simplified analytical version 
of the abrasion-ablation model which is valid in regions 
where only neutron evaporation plays a role, has been 
shown to reproduce the cross sections of many neutron- 
rich fragments very well (see, e.g., Refs. [Dl[l6|,[l9j]). The 
same is true for the Z = 50 — 53 data points in Fig. [8] 
The model tends to slightly underpredict, however, the 
lowest-Z data points, where EPAX 3 seems to follow the 
data a bit better. 
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FIG. 6. (color online) Isotope 
distributions from 1 A GeV 
208 Pb fragmentation in a na 'Cu 
(circles, Ref. and a 9 Be 

(squares, Ref. [16|) target. The 
Cu-target data were scaled by 
a factor of 0.77 according to 
Eq. (|3]). The full line represents 
the new EPAX 3 parametriza- 
tion, whereas the dashed line 
shows the old EPAX 2 version. 
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FIG. 7. (color online) Same 
as Fig. [6] but for the system 



3 Xe+ a Be hj 



EPAX 3 be- 



haves worse than EPAX 2 very 
close to the projectile but fits 
better to the measured data 
points for most of the very 
neutron-rich fragments. 
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FIG. 8. (color online) Cross sections of proton-loss channels 
for 1 A GeV 136 Xe + 9 Be (dots, Ref. [H) in comparison with 
EPAX 3 (full line), EPAX 2 (dashed curve), and with the 
physical "cold fragmentation" model COFRA [2^| (stars). 



C. Fragmentation of proton-rich projectiles 

Contrary to the neutron-rich fragments discussed in 
the previous subsection, cross sections for even extremely 
proton-rich fragments like 48 Ni or 100 Sn were already 
predicted with surprising accuracy by EPAX 2. Care 
was undertaken in the present study not to deteriorate 
the excellent predictive power of EPAX on the proton- 
rich side. Fig. [9] shows that this has been achieved to a 
large extent. The figure shows measured cross sections 
for 112 Sn+ 9 Be [15| . complemented by data for less ex- 
otic fragments from the system 112 Sn+ 112 Sn [26(, both 
data sets obtained at 1 A GeV. The Sn-target data were 
scaled by a factor of 0.65 (see Eq. ([3])) to match the Be- 
target data. Fig. [S] confirms that both data sets fit to- 
gether smoothly and that EPAX 3 describes the entire 
distributions rather well, very similar to EPAX 2. As a 
benchmark test, one can check the cross section of 100 Sn 
from 112 Sn+ 9 Be fragmentation: The agreement with the 
measured value of l.Sj^jj pb [ll| has deteriorated a bit, 
but not much (Version 3 predicts 7.9 pb, compared to 
6.6 pb for Version 2). Another benchmark is provided 
by the fragment 48 Ni produced from 58 Ni+ 58 Ni at 160 A 
MeV with a cross section of 100 ± 30 fb [28]. This is to 
be compared with an EPAX 3 prediction of 20 fb, EPAX 
2 predicted 57 fb. 



D. Range of validity of the EPAX formula 

As mentioned in the Introduction, the EPAX 
parametrization contains no bombarding-energy depen- 
dent terms, following the concept of "limiting fragmen- 
tation", which assumes that the cross sections for the 
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fragment mass A 

FIG. 9. (color online) EPAX 3 (full lines) and EPAX 2 
(dashed lines) predictions for production cross sections from 
112 Sn+ 9 Be, in comparison with data on a 9 Be (Ref. [T^ |. 
squares) and a 112 Sn (Ref. [2q |. circles) target measured at 
1 A GeV. The latter data points have been scaled by a factor 
of 0.65. 



formation of fragments rather close to the projectile are 
independent of bombarding energy if the latter is suffi- 
ciently above the Fermi energy in nuclei (« 40^4 MeV). 
This concept is difficult to prove since only few systems 
have been investigated where identical fragments were 
produced at several different bombarding energies. One 
system is Ni+ Be which was studied at 140 and 650 
A MeV [ID, [2(| • An example how well results from 
these two studies fit to each other was shown above in 
Fig. [T] other mass chains exhibit similar agreement. An- 
other study also points to the fact that around 140 A 
MeV energy independence is largely achieved: the re- 
action 76 Ge+ 9 Be at 132 A MeV studied by Tarasov et 
al. [HI . Most of their cross sections of very neutron-rich 
fragments in the range 23 < Z < 17 are quantitatively 
reproduced by EPAX 3. 

On the other hand, a comparison of two 86 Kr-induced 
data sets, one measured at 500 A MeV [23[ and the other 
at 64 A MeV [27j], do not yield identical results. EPAX 
fits to the latter yield systematically narrower width pa- 
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rameters, R. Moreover, the simple geometrical scaling 
law for two different targets, which seems to work at 
high energies (see Figs. l6l9| ). breaks down as visualized in 
Fig. 11 of Ref. [27). In agreement with Mocko et al. we 
conclude that an energy of 64 A MeV is below the range 
of validity of the EPAX formula. 

IV. SUMMARY 

Some modifications of the EPAX 2 parameters were 
shown to improve the predictive power of the EPAX for- 
mula in regions where too large cross sections were calcu- 
lated, in particular for very neutron-rich fragments from 
medium-mass and heavy projectiles. In addition to a 
modified parameter set, a "downscale factor" had to be 
introduced for neutron-rich projectiles that depends on 
the neutron-excess of the projectile and scales with the 
neutron-excess of the fragment. The resulting EPAX 3 
formula reproduces with rather good accuracy the bulk of 
the measured fragmentation cross sections as long as the 
incident projectile energy is larger than 130 - 140 A MeV. 
Systematic discrepancies were found for 86 Kr projectiles 



at 64 A MeV. The improvement of EPAX 3 compared 
to the previous version is confirmed by the much better 
agreement with the physical "cold-fragmentation" model 
COFRA, which was found to reproduce well measured 
cross sections for the most neutron-rich fragments from 
136 Xe and 208 Pb. 

The program code of the EPAX formula can be down- 
loaded from the GSI Document Server [291 ] . 
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